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The fully covariant formulation of transformation optics is used to find the configuration 
of a cloaking device operating in an expanding universe modelled by a Friedmann-Lemaitre- 
Robertson-Walker spacetime. This spacetime cloak is used as a platform for probing the covariant 
formulation of transformation optics, thereby rigorously enhancing the conceptual understanding 
of the theory. By studying the problem in both comoving and physical coordinates we explicitly 
demonstrate the preservation of general covariance of electrodynamics under the transformation 
optics procedure. This platform also enables a detailed study of the various transformations that 
arise in transformation optics. We define a corporeal transformation as the “transformation” of 
transformation optics, and distinguish it from coordinate and frame transformations. We find that 
corporeal transformations considered in the literature have generally been restricted to a subset of 
all possible corporeal transformations, providing a potential mechanism for increased functionality 
of transformation optics. 
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I. INTRODUCTION 

The prediction and subsequent laboratory demonstration of an electromagnetic cloaking device has created a fun¬ 
damentally new approach to the design of electromagnetic and optical device^ called “transformation optics” or 
“transformation electrodynamics” (which we shall generically abbreviate as TO)[ll-lj. The approach has been widely 
explored over the last decade both theoretically and experimentally, with studies of cloaks of various shapes, sizes, and 
operating frequencies; and applications beyond cloaking to myriad different optical systems including lenses, waveg¬ 
uides, antennas, and beam splitters. The transformation method has even been extended beyond electrodynamics to 
encompass acoustics, elastics, thermodynamics, and even quantum systems. Several extensive reviews are available 
which delve into the increasingly voluminous literature related to cloaking and the transformation method, including 
refs. [1,0. 

The TO design paradigm introduced by Pendry et. al. [ij was based on the covariance of Maxwell’s equations under 
coordinate transformations. In a nutshell, shifting from one coordinate system to another changes the coefficients in 
Maxwell’s equations, but it was found that these transformed coefficients could be absorbed by a redefinition of the 
permittivity and permeability within the original coordinates. The redefined material parameters essentially introduce 
a non-vacuum medium in place the of the coordinate transformation. 

Since a coordinate transformation is a diffeomorphism that leaves the system physically invariant, it is unclear from 
this nutshell description how TO actually works. Potential misunderstanding stems from the fact there are several 
different mathematical operations involved, all of which are usually referred to as some kind of “transformation”. Even 
more unfortunately, they are all likely to be referred to as “coordinate transformations,” and significant confusion 
exists within the field as to the nature and meaning of these different transformations. 

In ref. 0, the idea of passive vs. active transformations was discussed. A passive transformation is a true coordinate 
transformation that leaves the physical object (a vector, say) unchanged while, for example, rotating the coordinate 
system. The same vector has different component values relative to the original and rotated coordinate systems, but 
still represents the same object. An active transformation, on the other hand, actually changes the physical quantity 
while leaving the coordinates unchanged - rotating the actual vector rather than rotating the coordinates. 

In TO, the main transformation at work is one that results in the description of a physical object ~ a “transfor¬ 
mation medium” - whose parthenogenesis from the vacuum must be in some way related to the notion of an active 
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transformation rather than that of a passive coordinate transformation. The word corporeal means “having a material 
existence” or “consisting of a material object;” and in the spirit of this word we introduce the terminology corporeal 
transformation to mean a transformation that results in the description of a transformation medium. Thus the “trans¬ 
formation” of transformation optics is a corporeal transformation, rather than a true coordinate transformation. In 
addition to corporeal and coordinate transformations, we will find that frame transformations also play an important 
role in TO. 

In this paper we study cloaking in a curved spacetime, specifically in an expanding Friedmann-Lemaitre-Robertson- 
Walker (FLRW) spacetime. Friedmann-Lemaitre-Robertson-Walker has some physical relevance as the spacetime of 
modern cosmological models describing our expanding universe, although the actual cloak parameters for such a cloak 
are of limited practical value, and obtaining the FLRW spacetime cloak parameters is not the primary objective of 
this paper. Instead, the FLRW cloak serves as a vehicle for investigating the underlying framework of TO, and there 
are three main results of this investigation: 


1. We show that general covariance is carried across the corporeal transformation, meaning that an actual coor¬ 
dinate transformation may be applied either before or after the corporeal transformation without changing the 
physical nature of the transformation medium. Although one hopes this should be true by the construction of 
the general theory used here it is not entirely obvious and not guaranteed. In fact it is not possible to rigor¬ 
ously answer this question within the more widely used constructions of transformation optics that interpret 
the transformation as a coordinate transformation associated with a fictitious curved spacetime. 


2. The distinction between corporeal, coordinate, and frame transformations is clarified through a detailed study 
of the FLRW cloak. The FLRW cloak is the ideal platform for this work because it is the simplest cloak and 
spacetime for which the distinctions become readily apparent and most easily explicable. 


3. We find that the set of all corporeal transformations may be divided into two classes which we will call “coordinate 
corporeal transformations” and “frame corporeal transformations.” The frame corporeal transformations present 
a new avenue for further study in TO and may be an interesting and useful class of corporeal transformations. 


The most natural framework for TO is that of general relativity, where light propagates on a 4-dimensional, possibly 
curved manifold. The connection between TO and curved spacetime was first hinted at in ref. Q where the propagation 
of light through a fictitious curved spacetime was identified with the propagation of light through a medium residing 
in flat spacetime. But that approach was both restrictive in it’s scope and included an artificial and non-rigorous 
interpretation of a coordinate transformation of flat spacetime as a curved spacetime. Later it was shown how TO 
may be more rigorously formulated in such a way that the spacetime manifold under consideration is always the 
spacetime manifold where the device will operate, enabling the TO paradigm to be applied in curved spacetime 
environments Although for many years general relativity was seen as an esoteric subject with no real-world 
applications, it is well known that general relativistic corrections are crucial to the operation of global positioning 
systems, and so may not be ignorable for future near-Earth TO applications. Such corrections have been studied 
for TO in the Schwarzschild spacetime of a massive spherical object such as Earth Q. There is, therefore, merit in 
developing the most complete and general theory possible and in exploring its capabilities and limits. Eriedmann- 
Lemaitre-Robertson-Walker spacetime is a nice choice for this study because while it is relatively simple it exhibits 
rich physics. 

This paper is organized as follows: In Sec.|II]we review the completely covariant, 4-dimensional theory of electrody¬ 
namics and dielectric media on spacetime manifolds. This form of electrodynamics leads naturally to the formulation 
of TO described in Sec. mi In Sec. m the FLRW spacetime is described, the appropriate cloak transformations 
are defined, and the covariant theory of TO is employed to determine the cloak parameters. The resultant cloak 
is subjected to coordinate transformations and studied in the flat space of a local observer, demonstrating that the 
general covariance of electrodynamics is preserved under the operation of TO. In Sec. |V] we interpret the results of 
Sec.HV] and suggest how corporeal transformations considered in the literature to date have generally been restricted 
to a subset of the possible corporeal transformations, providing a potential mechanism for increased functionality of 
the TO paradigm. A very detailed analysis of light propagation in FLRW spacetime is provided in Appendix]^ while 
detailed discussions of coordinate transformations and frame transformations are relegated to Appendices [B1 [Cl and 
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II. COVARIANT ELECTRODYNAMICS 

The standard vectorial representation of Maxwell’s equations 
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makes the assumption of flat, Minkowski spacetime. But it has been shown that transformation electrodynamics can 
only be fully understood in the context of a theory that can adequately distinguish spacetime from medium ili3, 
so it is advantageous to work with a formulation of electrodynamics where this distinction can be easily made, even 
if the spacetime is flat. Here we review the salient features of a formulation that is manifestly covariant, containing 
only tensorial objects and the spacetime metric [^[^. 

We let spacetime consist of a manifold M along with a symmetric metric g that describes the shape of the spacetime. 
Given a curve on the spacetime manifold, the squared infinitesimal distance between two points on the curve is given 
by the line element 


ds^ = g^udx^dx’' ( 2 ) 

where is the matrix representation of g. We adopt the sign convention (—h H—h) and we use units where the 
speed of light c = 1. With these conventions the metric of Minkowski spacetime in Cartesian coordinates is 
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Although the scalar and vector potentials are often treated as non-physical mathematical tools, in light of their 
role in quantization of the field we follow Maxwell himself in maintaining their central role, although now combined 
into the 1-form potential A = where = (/i. A) in Minkowski spacetime. The field strength tensor F = F^i, is 
obtained by taking the exterior derivative of A 

F = dA, = d^A, - d,A^ (4) 


This is the 4-dimensional generalization of obtaining E and B as derivatives of the potentials, so it follows that F 
contains the field components. In the locally flat space of an observer, the components of F form an antisymmetric 
matrix 
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On the other hand, we let the components of the electric displacement D and magnetic field F[ be contained in the 
excitation tensor G. In a locally flat space, the components of G form an antisymmetric matrix 
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As they stand. Maxwell’s Equations (Hj) are incomplete and require the specification of an independent relationship 
between the fields before solutions can be found, and the same holds true when the fields are expressed in tensorial 
form. In the Minkowski vacuum, that relationship reduces to equating components Fla = Ba and Da = Ea- This can 
be made more precise by employing the Hodge dual *, a well-defined and natural operation on pseudo-Riemannian 
manifolds endowed with a metric, such that in vacuum 


G = *F 


( 7 ) 
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or, in component form 


G;.. = (*F)^, = ^ Vl^ . (8) 

Consider Earth’s oceans for a moment. A dense object dropped into the ocean will sink, experiencing the inexorable 
pull of gravity regardless of how much water surrounds it. But gravity is nothing more than a word used to describe the 
curvature of spacetime, and from the ocean example it is clear that the spacetime surrounds and permeates everything 
that exists; it cannot be removed nor ignored. Realizing that dielectric media exists within, and is permeated by, 
the spacetime, we seek a formulation of Maxwell’s equations that clearly distinguishes any spacetime or metric 
contributions from the medium contributions. It has been shown [^, Q that for linear media such a formulation can 
be obtained through a simple extension of the vacuum constitutive relation Eq. © to 


G = *xF (9) 

where the tensor x encodes parameters equivalent to permittivity, permeability, and magnetoelectric couplings. In 
component form the constitutive relation is 

( 10 ) 

The most general x has 36 independent components, and we require that i® independently antisymmetric on 

each pair of indices, a •<->■ /? and a p. Additional symmetry conditions may be imposed based on thermodynamic 
arguments, in particular that Xapap symmetric under the pair exchange a/? -O- ap. The general linear constitutive 


relation Eq. Q can be made applicable to the vacuum, thereby superseding Eq. by requiring 

XvacF = F, (II) 

which is sufficient to specify Xvac uniquely. Expanding Eq. (nni, the coefficients may be collected into the form 

Da = {e^),^Eb + {V)yb, = + ( 12 ) 

or in the more traditional form 

Da = eJ>Eh + Ba = pJ’Bb + (13) 

The 3x3 matrices e'’, etc. of the ‘covariant representation’ Eqs. (fT^ may be connected to the 3x3 matrices p, 
e, etc. of the ‘traditional representation’ Eqs. (HSl) by 

= (14) 

Finally, the covariant, 4-dimensional formulation of Maxwell’s equations may be summarized as 

dF = 0, dG = J, G = *xF, (15) 

where J is the free charge and current source, which in what follows we take to be zero. 


III. TRANSFORMATION OPTICS 

A good first step at generalizing transformation optics was to employ the Plebanski equations to identify a fictitious 
spacetime with a dielectric material in flat spacetime One disadvantage of the Plebanski based approach 

is that the final transformation medium is restricted to applications in vacuum Minkowski spacetimes. Since the 
Plebanski equations are furthermore not strictly covariant this approach fails to satisfy the desire for a completely 
general, manifestly covariant theory of TO. A strictly covariant theory is necessary for a complete and rigorous 
foundation, and to fully understand the role of transformations as studied here. Several authors contributed to the 
development of such a covariant theory (l2l - [l5l| . but here we review the approach initially presented in refs, 
as modified slightly in ref. 

The idea of TO can be formalized as follows. Begin with an initial configuration consisting of a spacetime manifold 
M and metric g, with associated Hodge dual *, and initial material distribution x (enabling TO within a prior 
non-vacuum medium 0). Let F and G be the electromagnetic field tensors on M, obeying Maxwell’s equations and 
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FIG. 1. Manifold M is mapped under T to the image M C M. T may be defined as a cloak transformation such that M 
contains a ’hole’ relative to M. However, the fields are transformed by the pullback of T = T~^. 


the constitutive relation G = Let T : M ^ M C M he a map of M to the image M, as in Fig. [TJ The image 
M may not cover all of M; for example, M may contain a hole relative to M, leading to an invisibility cloak. 

Given a map between two manifolds ip : M ^ N there exists an associated pushforward map that maps vector 
fields in the tangent space of M to the tangent space of N, and a pullback map ip* that maps covector fields from 
the cotangent space of N to the cotangent space of M. For differential forms like F the associated map of interest is 
the pullback. Since the pullback map ip* : T*N —)• T*M is opposite the direction oi ip : M ^ N itself, it follows that 
in TO we cannot use the pullback of T : M —M to map F from the cotangent space of M to the cotangent space of 
M. Instead we need the pullback T* of a map T '■ M M. The pullback map T* takes the form of the Jacobian 
matrix of T. 

The map T has taken center stage in TO because it provides an intuitive visualization of the desired propagation 
of light through the transformation medium and an easy method for picking a corporeal transformation based on a 
mapping of the trajectories of light - curves on the manifold ~ to new trajectories. In this picture an invisibility cloak 
can be imagined as mapping the space to a new space with a hole in it, but when calculating the transformation 
medium we actually require the inverse of T; in other words we may take 7” = T~^. In fact, TO may be pursued 
without recourse to the pullback T* at all. Instead, a map P : T*M —> T*M between cotangent bundles may be 
specified independently of a map between manifolds, as further discussed in the following sections. 

The essential feature is that the fields (F, G) are mapped to a new set of fields ^F, G^ on M, but for these new 

fields to be valid solutions of Maxwell’s equations we require a new ‘transformation’ medium x- Furthermore, since 
TO can not change the spacetime, only the medium, M C M inherits the metric from M. From the constitutive 
relations on M and M we find 


G = T*iG)=r* ix*F) = x*T*F 


which can be solved for the transformation medium @, i [13,13 
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In Eq. m, A is the Jacobian matrix of T, A ^ is the matrix inverse of A. The initial material tensor x ^-nd the 
first * are evaluated at T{x), while everything else is evaluated at a; G M. 


IV. CLOAKING IN FLRW SPACETIME 
A. ELRW spacetime 

In 1929, Edwin Hubble discovered that not only were all distant galaxies receding from us but that the recessional 
speed of more distant galaxies was faster than the recessional speed of less distant galaxies. This phsenomenon, 
now referred to as Hubble’s Law, appears to be universal and holds true for any choice of origin - on large scales, 
everything in the universe is moving away from everything else. The interpretation of this Hubble flow is that the 
universe is uniformly expanding, with all the galaxies we see going along for the ride. Not all large-scale motion is 
due solely to the universal expansion, but as reviewed in Appendix [^ this extra peculiar velocity is damped out as 
the universe expands, so that eventually all motion will asymptote to the Hubble flow. Such an expanding universe 
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is well-described by a Friedmann-Lemaitre-Robertson-Walker (FLRW) spacetime, which is the general relativistic 
solution of a spacetime filled with a perfect fluid. 

The FLRW spacetime is typically described in “comoving” spherical coordinates {t,rc,d,ip), with the line element 


and corresponding metric 


= —dt^ + a(t)^ {dri + rldfl^) 
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where a{t) is a scale factor that governs the size of a spatial slice of the manifold. From the line element one can see 
that if a{t) is either increasing or decreasing with time then the spatial universe is either expanding or contracting. 
For definiteness, in what follows we assume the universe is expanding. Comoving coordinates are somewhat peculiar 
in that they are fixed to the spacetime. Since the spatial universe is expanding, it follows that the physical separation 
between two points of fixed comoving coordinates is actually increasing with the expansion, similar to the way that 
the distance between two nearby dots drawn on the surface of a balloon will increase as the balloon is inflated. An 
object at rest at radius in comoving coordinates has an increasing physical distance d = a(t)rc. Objects at rest in 
comoving coordinates are said to be comoving, and a comoving object has 4-velocity 


ui^ = (1,0,0,0). 


( 20 ) 


On the other hand, an observer at the origin would interpret the recession of objects due to universal expansion 
as some non-zero velocity of the object. To reflect this observer experience, it is useful to describe FLRW spacetime 
in “physical” coordinates {t,rp,6,(p). Physical coordinates are related to comoving coordinates via the coordinate 
transformation 


Tp = a{t)rc, 


( 21 ) 


and in these coordinates the line element becomes 

a(tf 


ds^ = -\l- 


- 2^^^rpdrpdt + dr^ + rpdfl'^ 
where an overdot denotes the derivative with respect to t, with corresponding metric 


( 22 ) 
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In this and what follows, primed indices refer to physical coordinates while unprimed indices refer to comoving 
coordinates. In physical coordinates, a comoving object has 4-velocity 


n^' = (l,-rp,0,0), (24) 

a 

with the characteristic feature that the radial velocity depends on the distance, in agreement with Hubble’s law and 
providing compelling evidence that the universe is well-modelled by an expanding FLRW spacetime. On the other 
hand, an object at rest in physical coordinates has normalized 4-velocity 


u 



(25) 


The comoving and physical pictures are two equally valid coordinate descriptions of the same spacetime, with a 
well-defined coordinate transformation to switch between the two. Despite the fact that the different coordinates 
describe the same physics, the geometric optics description of light propagation appears quite different for an observer 
using comoving vs. physical coordinates. By solving the geodesic equation for null geodesics in each case, as outlined 









7 


in Appendix parametric equations of null trajectories are found, with the result plotted in Fig. In the figure, 
the affine parameter is increasing from left to right, and in Fig. |2(b)| the initial conditions specify a congruence of 
initially “parallel” geodesics. It is clear from the figure that the notion of parallel does not have its usual meaning in 
this space of uniform curvature (constant curvature over a constant-time spatial hypersurface). 

Note that by making a change of time coordinate, Eq. can be rewritten as 

= a(r)^ (—-I- dr^ + , (26) 


which is just a conformally scaled version of Minkowski space. Since conformal maps preserve angles the null geodesics 
of Minkowski spacetime are conformally mapped to null geodesics in FLRW, thereby preserving the path (although the 
propagation along the path is different). Thus it is no surprise that Fig. 2(a) looks just like Minkowski spacetime, but 


we will see that looks can be deceiving when it comes to defining a cloak transformation in the comoving coordinates 
of FLRW spacetime. 



(a) Comoving coordinates (b) Physical coordinates 


FIG. 2. Null curves of FLRW spacetime plotted in (a) comoving coordinates, and (b) physical coordinates. Despite the non¬ 
zero curvature of FLRW spacetime its conformal flatness is apparent in comoving coordinates, where initially parallel geodesics 
remain parallel as in Minkowski spacetime. The expansion of an initially parallel congruence of null geodesics is made clear in 
physical coordinates. 


B. Cloak transformations in FLRW 

Given the multiple coordinate choices of FLRW spacetime and the different propagation pictures provided in each, 
the next step is to choose an appropriate corporeal transformation. Consider first the Pendry et. al. transformation 
[I| in comoving coordinates 
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Bear in mind that this is part of the map T : M —>■ M C M; the complete specihcation of T is 
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Also recall that to calculate the transformation medium x with Eq. (11711 we actually require the inverse T = T~^ 
rather than T itself. 

This transformation maps the origin to the cloak’s inner surface at fixed radius ?"c(l), but in comoving coordinates 
the fixed radius ^^(l) defines the surface of a physically expanding sphere. To a local observer, such a comoving cloak 
would appear to be expanding with the expansion of the universe. Although one could, with all validity in the context 
of transformation optics, use this corporeal transformation and calculate the corresponding transformation medium, 
a physically expanding cloak is not what we had in mind. 

Instead, the picture we have in mind is of a cloak of fixed physical size relative to a real, physical observer. Every 
element of a cloak of fixed physical size with inner radius i?i and outer radius i ?2 would actually have some non-zero 
velocity in comoving coordinates. This means that the cloak medium would not be comoving, being instead held 
together by whatever electrostatic forces are present to overcome the universal expansion and keep the cloak from 



















































breaking apart. In other words, a cloak of fixed physical size would actually appear to be shrinking if described in 
comoving coordinates. 

On the other hand, in physical coordinates a cloak of fixed physical size would not appear unordinary to a real, 
physical observer, and there the transformation 

Tp —>■ i?i H-—- Vp (29) 

rt2 

would be appropriate. It is not immediately obvious whether the Pendry transformation will work in this setting of 
universal expansion, but the action of a cloak derived from this corporeal transformation in physical coordinates is 
depicted in Fig. O and it may be observed that the downstream ray trajectories are perfectly preserved. 


yp 



FIG. 3. Equatorial slice of a cloak in physical coordinates of FLRW spacetime. The universal expansion of a set of initially 
parallel null curves is accommodated by the cloak transformation. 


Using the coordinate transformation Eq. m, the corporeal cloak transformation Eq. may be transformed to 
comoving coordinates, where the corporeal transformation becomes 


Ri R2 — Ri 


(30) 


which for increasing a{t) does indeed represent a shrinking cloak transformation, as predicted for comoving coordinates. 

As will be shown, Eqs. (EH) and (l30l) are corporeal transformations for the same cloaking device described by two 
different coordinates systems. Recalling that it is the inverse transformation that appears in Eq. we find the 
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in physical coordinates. 

Using the cloak transformations Eqs. (|30)) and (l29ll . their associated Jacobian matrices Eqs. EID and (l32|) - and 
the metrics Eqs. m and (|23l) . the tensorial components of the transformation medium x may be calculated in each 
coordinate system through the application of Eq. dUD (we assume the initial space is vacuum, x = Xvac)- We write 
the components as and i^p)^ where again unprimed and primed indices refer to comoving and 

physical coordinates, respectively. 
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C. Local Frame Results 


However, the coordinate components Xa'p'^ '' (^p) of somewhat limited value since they refer 

to coordinate basis fields of the global coordinates, whereas engineers tasked with building the cloaking device would 
prefer to work in their locally flat frame. This is because a local observer measures electric and magnetic fields relative 
to her local frame, and the permeability and permittivity, as derived from x in 0 , only have their traditional meaning 
relative to that local frame. 

One of the defining features of a classically curved spacetime is that if you look at only a small portion of the 
manifold it appears to be Minkowskian, in the same way that the curvature of the Earth is ignorable on the local 
scale. By saying that the manifold “appears locally Minkowskian” we mean that it is possible, at each point, to define 
a basis {ba} of the tangent space and related basis {0^} of the cotangent space such that the line element relative to 
this basis is 

ds^ = (33) 


where rjAB is the Minkowski metric. This is similar to the way the line element takes on different forms in different 
coordinates, such as Eqs. (HU), (USD, and (E51) . However, whereas the tangent and cotangent space bases in Eqs. o, 
(1^ . and (E51) follow from the choice of coordinates, the so called “coordinate bases,” here the tangent and cotangent 
bases are changed independently of the coordinates, as further described in Appendices iBl and ICl The main point is 
that a basis may be chosen such that relative to the chosen basis the metric becomes Minkowskian even though the 
spacetime is fundamentally curved, which makes this the natural choice for observers working relative to their locally 
flat environment. 

We have seen that using the corporeal transformation Eqs. dSO]) and (1291) results in transformation media 
and Xa’B’^ ^ i^p)- described in Appendix [Cl the components of x in a local frame of comoving coordinates are 


XAB^^i^c) = eA°'e/e^^e^^Xc,p^''i.Xc) 


while the components of x in a local frame of physical coordinates are 
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It still remains to find the frame transformation matrices e^“ and e^,“ . It may be seen from Eqs. (1551) and (IC9|) that 
the frame transformation matrix can be determined from the condition 


^A ^B dap = VAB 


(36) 


where 
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is the Minkowski metric in spherical comoving coordinates. An obviously similar condition holds for e^/“ in spherical 
physical coordinates. In additional to this condition, we demand that the 4-velocity of an observer at rest in the 
coordinate basis, e.g. Eqs. (gni) and (Ha, is frame transformed to the four velocity of an observer at rest in the local 
frame basis. This is equivalent to setting Cq = u^. One finds 
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in comoving and physical coordinates, respectively, with ^ and found from the inverse transpose of these 

matrices. The permeability and permittivity relative to the local observer may now be calculated as follows: The 
corporeal transformations in each coordinate systems are given by Eqs. (1301) and (1291) : the Jacobian matrices Eqs. dUD 
and (15^ of the corporeal transformations are used in Eq. m to calculate % relative to the coordinate basis; the frame 
transformation matrices Eqs. (l38l) and (l3^ are used in Eq. (1351) to calculate x relative to the locally flat frame; the 
coefficients of x are collected into three dimensional matrices through the identifications of Eq. (TTO)) (the component 
identifications may also be found in ref. 0), which are finally put into the traditional representation through Eqs. 
(HI. The result is 
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in a local frame of comoving coordinates, and 
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(R2-RlYa'^+Rl{rp-Ri)'^a'^ ) 


(40a) 


(40b) 


(41a) 


7 = 


-'‘7 = 


0 


0 

0 

R]_a (R2+Ri^p-^R2rp')a^+Rlrp(rp-Ri)d^ 
~ iR2-Ri)'^a‘^+Rl{rp-Riy^d‘^ 


Rid {Rl+Rirp-2R2rp'ja'^+Rlrp{rp-Ri)d.'^ 
iR2-Ri)'^a^+Rl{rp-Ri)'^d.'^ 

0 


(41b) 


in the local frame of physical coordinates. It may be readily seen that setting a{t) = 1 returns us to Minkowski 
spacetime where the results for both comoving and physical coordinates reduce to the ordinary flat spacetime Pendry 
cloak. 


D. Coordinate Transformations 


Finally, we can now show that the TO procedure and the resulting transformation medium is generally covariant, 
meaning that knowing the coordinate transformation to switch between comoving and physical coordinates, we can 
relate the coordinate expressions of x by (see Appendix [B]) 


■V ^ It f f f ^ -v ^ 

Xa'/3' y^p) — ^ a'^ P'Xa0 




^C{xp 


(42) 


where is the Jacobian matrix of the coordinate transformation Xp = C{xc) and is the Jacobian matrix of 
the inverse coordinate transformation Xc = C~^{xp). As applied to a tensor, the coordinate transformation is enacted 
by contracting the Jacobian matrices of the coordinate transformation and inverse coordinate transformation with 
the tensor. However, the Jacobian matrices of the coordinate transformation, by construction, refer to the coordinate 
basis of the tangent bundle, hence by themselves can only transform from the coordinate basis representation in one 
coordinate system to the coordinate basis representation in the other coordinate system. Instead, what we would 
really like is to transform from the local frame of comoving coordinates to the local frame of physical coordinates. As 
shown in Appendix [D1 the coordinate and frame transformations may be combined such that the components of x in 
one local frame may be directly transformed to the other local frame 


XA'B'^ ^ (2^p) — jj XAB^^i^c) 


■T[xp) 


(43) 
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with 


A _ T OL ^ < 

A' ^ oc^Oi'^A' 


(44) 


and 

M^'c = ec'^J^'e^V- (45) 

It is now possible to verify that Eqs. (14(111 correctly transform to Eqs. (l4T1l under the combined action of the coordinate 
and frame transformations embodied in Eq. (1431) . 

There is a widespread misconception within the field of TO that the formulation and results of transformation optics 
are intimately tied to a choice of coordinates and a subsequent coordinate transformation. Here we have explicitly 
demonstrated that this is not the case. A corporeal transformation given in one coordinate system may be given a 
different description in another coordinate system, and the resulting transformation media are physically equivalent 
and are mutually and directly related via a well-defined coordinate transformation. 

Although one hopes this to be true by the construction of Eq. dm), it turns out that it is not an entirely trivial 
or obvious result. The complete derivation of Eq. dm requires a more rigorous specification of the metric on each 
manifold Q; in particular, we may generalize the picture of Fig. [T] by mapping between different spacetime manifolds. 
This generalization has been used to describe dielectric analog spacetime models that mimic the propagation of light 
through a vacuum spacetime by the propagation of light through a suitably chosen dielectric medium, and forms 
the basis for a rigorous derivation of Plebanski’s equations [13. However, in that case the coordinates of the target 
manifold are physically distinct from those of the domain, thereby breaking strict covariance. Thus, for example, 
the FLRW spacetime described in either comoving or physical coordinates may be mapped into an analog medium 
described by spherical coordinates of Minkowski spacetime, but these will be physically inequivalent representations of 
the spacetime that cannot be related through a simple coordinate transformation. The reason covariance is maintained 
across the corporeal transformation is because in TO we have insisted on the preservation of the underlying spacetime. 
Thus it would seem that formulations of transformation optics based on the Plebanski equations, or equivalently on 
the analog model, are incapable of rigorously demonstrating the general covariance of the transformation media. The 
lack of strict covariance in analog models will be discussed in more detail elsewhere. 


V. DISCUSSION AND CONCLUSIONS 

We have studied cloaking in an expanding Friedmann-Lemaitre-Robertson-Walker spacetime, which is a nice choice 
because while it is a relatively simple spacetime it is rich in physics and serves as the model for our expanding universe. 
The FLRW cloak makes an ideal platform from which to examine and clarify the role of several different kinds of 
transformations that may appear in TO. In particular it was shown how medium-producing TO transformations, here 
referred to as corporeal transformations, are distinct from coordinate transformations and frame transformations, and 
it was shown that all three types of transformations play distinct but important roles in TO. The corporeal transfor¬ 
mation is responsible for what manifests itself as a transformation medium. Coordinate and frame transformations, 
on the other hand, simply provide different descriptions of the same physical phaenomenon and cannot be used to 
generate a new medium. A corporeal transformation may be referred to different coordinate systems, but the physics 
carried by the corporeal transformation remains invariant under the coordinate transformation, resulting in different 
coordinate descriptions of the same transformation medium. In other words the transformation medium obtained 
from a corporeal transformation in one coordinate system is physically equivalent to the transformation medium 
obtained when the corporeal transformation is expressed in a different coordinate system, and the two descriptions of 
the transformation medium are covariantly related by the coordinate transformation. 

The FLRW cloak provides a good example of how deciding on the appropriate cloaking transformation depends on 
the point of view of the observer and their choice of coordinates, but that ultimately the same cloak transformation may 
be referred to different coordinate systems. The desirable scenario of a cloak of fixed physical size was obtained from 
the usual cloak transformation in physical coordinates, but when referred to comoving coordinates the transformation 
actually actually describes a radially shrinking cloak. The cloak transformations in different coordinates are related by 
the coordinate transformation, and it turns out that this relationship is preserved across the corporeal map such that 
the transformation media obtained from either coordinate system are physically equivalent and also related by the 
coordinate transformation. Although this coordinate invariance is clearly desirable it is not obvious, and in fact hinges 
on the fact that the corporeal transformation acts only on the fields, while leaving the manifold untouched. This 
analysis clearly shows that the general covariance of electrodynamics is preserved through the corporeal transformation 
into the transformation medium, but this need not have been guaranteed. In fact, the closely related formulation of 
analog spacetimes [13 from which the Plebanski equations may be derived and were used in earlier formulations of 
TO [7|, breaks covariance under the analog map. 
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A coordinate transformation acts on a tensor via the differential of the coordinate map. Practically, this differential 
manifests itself as the Jacobian matrix of the transformation, with one index referred to each coordinate system, 
and is implemented via matrix multiplication; one Jacobian matrix for each tensor index. The implementation of a 
frame transformation is also basically matrix multiplication on each tensor index, but while one index of the frame 
transformation matrix refers to the coordinates, the other index refers only to an arbitrary basis vector field that 
is not, in general, associated to any coordinate system. As described in Appendix [Cl a frame field assigns a set of 
basis vectors at every point, so a frame field is nothing more than a set of four vector fields that may either be chosen 
freely or chosen to satisfy some conditions, and any vector field can be decomposed relative to this chosen set of basis 
vector fields. Here the frame field is chosen such that the associated metric is Minkowskian. The arbitrary index on 
a frame transformation matrix refers to these basis vector fields. 

If the frame transformation is smooth then one may find integral curves of the frame basis fields. If, in addition, the 
frame fields {ca} commute, [e^, es] = 0, then the parametrized integral curves may be identified as a set of coordinate 
functions {y^} such that ba = and = dy^ act as mutually dual bases for the tangent and cotangent bundles. 
What this shows is that for commuting frame basis vector fields, the arbitrary index on the frame transformation 
matrix actually refers to the coordinates corresponding to the integral curves of the fields, and in this case the frame 
transformation is nothing more than a coordinate transformation. This provides a refined perspective of the Jacobian 
matrix, where instead of considering the indices as referring to coordinates^ they should be thought of as referring to 
the vector fields whose integral curves correspond to the coordinates, i.e. to the coordinate vector fields. With this 
refinement, the Jacobian of a coordinate transformation is actually just a special type of frame transformation matrix. 

Transformation optics was initially conceived as a coordinate transformation but we have shown that corporeal 
transformations are distinct from coordinate transformations. Nevertheless, there is a clear analogy in that each can 
be identified with a set of functions on the manifold that act upon tensors via the Jacobian matrix of those functions. 
On the other hand, we have seen how coordinate transformations of tensors are a special case of frame transformations; 
in other words, if {IH} is the set of all possible frame transformation tensors, then the coordinate transformations 
are the set of all C S 91 such that C = dT for some function T. It follows from the analogy with coordinate 
transformations that the space of possible corporeal transformations is actually much larger than what has generally 
been considered in the literature. In particular, in analogy with a frame transformation, a corporeal transformation 
may more generally be given through the direct specification of the transformation tensor. Letting {©} be the set 
of all possible corporeal transformation tensors, then the set of all U G 6 such that D = dT for some function T 
could reasonably be called “coordinate corporeal transformations.” Although coordinate transformations are actually 
a subset of frame transformations, the terminology “frame transformation” usually refers to frame transformations 
that do not belong to the subset of coordinate transformations. By analogy, we may refer to the set of corporeal 
transformations that are not coordinate corporeal transformations as “frame corporeal transformations.” 

With few exceptions the multifaceted “transformation physics” program has restricted itself to coordinate 

corporeal transformations and very little seems to have been done with more general classes of transformations. Frame 
corporeal transformations therefore present a potential mechanism for increased functionality of the TO paradigm. 
A detailed analysis of these frame corporeal transformations in the context of TO presents a new avenue for further 
study that lies beyond the scope of the present paper, but is the subject of ongoing investigation. 


Appendix A: Affine parametrization in spacetime 


By applying the map T (see, e.g. Fig. [T]), a parametrized curve x^(r) C M will be mapped to the parametrized 
curve x^(t) C M, which is the path followed by light traversing the transformation medium. Such direct visualization 
is one of the alluring features of transformation optics, in that the corporeal transformation can be elucidated from 
a deformation of an initial set of trajectories of light into a desired set of trajectories. However, to take advantage 
of this procedure it is necessary to know the initial set of trajectories. Since the parametrization of null geodesics 
on spacetime manifolds may be unfamiliar, here we explain how to obtain the parametrization of null geodesics in 
Friedmann-Lemaitre-Robertson-Walker spacetime. We do this in some detail because a complete derivation is not 
readily found in the literature. 

We parameterize a curve with affine parameter r, along with its velocity and acceleration as 


x^{t) = [x^( t), x^( t), x^( t), x^( t)) (Ala) 

v^{t) = —— = x^ = (x^(t),x^(t),x'^(t),x^(t)) (A-lb) 

= x>" = (ai° (r), T (r), x^ (r), x^ (r)) (Ale) 



where 
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,n, X dx^ , d?x^ , . 

X (t) = x°{t) = etc. (A2) 

It should be mentioned that although the point x^ = {x^, x^, x'^, x^) is often thought of as a “vector,” such a position 
vector is only a point on the manifold and is not an element of the tangent space, as is. In this setting, x^ is not 
a vector and does not transform like a vector. Thus a curve x^(t) transforms with the map T itself, while the vector 
field transforms via the pushforward of T (assuming it exists). 


1. Null Geodesics in Minkowski Spacetime 

Before considering null geodesics in Friedmann Robertson Walker spacetimes, it is worthwhile to review the explicit 
affine parameterization of null geodesics in Minkowski spacetime, first in Cartesian coordinates and then in spherical 
coordinates. 


a. Cartesian coordinates 


In Cartesian coordinates the line element for Minkowski spacetime is 

ds^ = —c^dt^ + dx^ + dy^ + dz^. (A3) 

For a null geodesic we require 

gapv°‘v^ = -cH^ +x^ +f + = Q. (A4) 

The geodesic equation is 

i^^+TlpX^xP, (A5) 

but in Cartesian coordinates F((^ = 0 so we are left with the set of equations x^ = 0. These integrate immediately to 

(t(T), x{t), yir), z{t)) = (wt, Ux, Uy, Uz)t + (to, xo, yo, ^o) (A6) 

with the requirement that u^ + Uy + = u^c^, and we generally set c = 1. 


b. Spherical polar coordinates 


Integrating the geodesic equations in Cartesian coordinates was immediate, but things get messy in spherical polar 
coordinates. The Minkowskian line element is now 

ds^ = -c^dt^ + dr^ + r^dO^ + sin^ 9 d(p^ (A7) 

and the null condition is 

ga/3V°‘v^ = -c^P + P + sin^ 0(/3^) = 0. (A8) 

The geodesic equations are now 

t = 0 (A9a) 

f — riiP' + sin^ = 0 (A9b) 


9 + 2-9 — cos0sin0(^^ = 0 
r 


(A9c) 
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(p + 2-ip + 2 cot 99ip = Q. (A9d) 

r 

The t equation of course integrates immediately to the same result as previously obtained. This should be expected 
since transforming between Cartesian and spherical coordinates leaves t unchanged. Although the ip equation can 
readily be integrated to 




UJo 


/ rp sin 6>o 
y r sin 9 


2 


(AlO) 


we quickly run into difficulties finding an analytic solution to the 9 equation. Things can be simplified by restricting 
motion to the equatorial plane, in which case 9 = 9q = tt j2. With this restriction, 


p = UJQ 



and the radial equation becomes 


which, with 


and the null condition 


f — 


2 4 
^3 


= 0 


■f(r = 0) = Ur, 


2 I 2 2 2 2 

Ur+UJo^O = UtC . 


can be directly integrated to 


(All) 


(A12) 


(A13) 


(A14) 


r(r) 


UtC 




1/2 


(A15) 


It is easy to verify that for purely radial geodesics, in which case wp = 0 and Ur = c, r(T) = rp ± cr as expected. 
Knowing r(r), Eq. (lAllI) can now be integrated to find the parametrization of ip 


ipir) = ± tan 


-1 


rpM, 




roUJQ 


=F tan 


-1 


T’OWp 


■ /?o- 


(A16) 


Choosing an initial point (tp, rp, 7r/2, (/3p) and initial velocity with respect to r (m*, m,., 0, wp) we then have an 
explicit expression for the curve. Furthermore, one can readily convert this to Cartesian coordinates with x(t) = 
r{T) cos {ip (t)), and ?/(t) = r(T) sin((/j(r)) if, for example, one wanted to make a parametric plot with something like 
Mathematica. These become 


X = {ur cos ipQ — rpwp sin (/?p)r + rp cos (/?p 
y = {ur sin t/jp + rpwp cos /5p )t + rp sin ipQ. 


(A17) 


Note that if the chosen initial point on the trajectory coincides with the point of closest approach, then xq = b 
is just the impact parameter, Ur = 0 because at the point of closest approach the velocity is purely tangential, 
and ipo = ipin + 7r/2 where ipin is the angle at which the trajectory crosses the x-axis (special care is needed for 
trajectories with initial point at the origin). If one is unconcerned with following the wavefront and only cares about 
the geometrical path, then rescaling r xquiqt allows us to write 


X = T cos ip in — b sin ipin 
y = TSimpin +bC0Sipin- 


(A18) 
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2. Null geodesics in FLRW: comoving Cartesian coordinates 

a. Null geodesics 

In comoving Cartesian coordinates the FLRW metric is 

ds'^ = —c^dt^ + a^{t) {dx^ + dy^ + dz'^) . (A19) 

We parameterize a curve with affine parameter r, along with it’s velocity and acceleration, in the same manner as 
Eqs. (|Alall - (|Alcp . It should be borne in mind that the scale factor a{t) appearing in the metric, which is a scalar 
function, is different from the acceleration vector a^. Note that with this parameterization, the scale factor is also 
parametrized by r through its dependence on t, in other words a = a(t(r)). We will suppress the explicit dependence 
on t or r except where necessary for clarity. Dot notation will denote derivatives with respect to r and primes will 
denote derivatives with respect to t 


For a null geodesic we require 


a' = a' {£) 


da 
dt ’ 


da . ■ 
a = — = at. 
dr 


gapv v>- 


= -cH"^ + a?{x^ + Z^) = 0 . 


(A20) 


(A21) 


The geodesic equation yields 



(A22a) 


2a' . 

X H- ix = 0. (A22b) 

a 

The equations for y and z are identical to that of a;, thus it is sufficient to consider only the equations for t and x. 
Using the null condition Eq. (|A2ip and a't = d, the t and x equations become 


d ■ 

t T —t = 0, 
a 


(A23a) 


X + 2-x = 0 
a 


with solutions 


t{T) = utao 



a 


(A23b) 


(A24a) 


r;(T) = u^al j ^ 


—dr. 


(A24b) 


and similarly for y and z. The coefficients {ut,Ux,Uy,Uz) give the direction of propagation of the geodesic at r = 0, 
i.e. = {ut,Ux,Uy,Uz), and satisfy the null condition 

— c A Uq(u 2 ; a Uy A = 0. (A25) 


b. Conformal flatness 


It is well known that FLRW is conformally fiat, which means that FLRW is equivalent to a conformal rescaling of 
Minkowski space. In other words, transforming to “conformal time” 


drj 



(A26) 
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the metric becomes 


ds'^ = a{t) (^—<?drf + dx^ + dy^ + dz^) . (A27) 

Since the null character of a tangent vector (and hence the null character of a geodesic) is invariant to conformal 
transformations, null geodesics of Minkowski spacetime are conformally mapped to null geodesics of FLRW spacetime. 
In Minkowski spacetime, a defining characteristic of null geodesics is that Ax = c At for light propagating along the 
x-axis, and this relation between intervals should also be reflected in FLRW if we switch to conformal time. Indeed, 
transforming Eq. (IA23al) to conformal time, we find the equation 

I;+ 2^77 = 0 (A28) 

with corresponding solution 

77 ( t ) = UrO-l J (A29) 

Thus the intervals Ax = c Ay clearly show the conformal flatness of FLRW. 


3. Evaluating the scale factor 


To proceed further requires knowledge about the form of the scale factor function a(r). The scale factor cannot 
be specified arbitrarily, but instead must be obtained from Einstein’s equations, It is known that the 

FLRW metric is the solution for a spacetime filled with a perfect fluid. The stress tensor for a perfect fluid with 
4-velocity is given by 

V 

= pUfj.u^ + + U^u^) (A30) 

where p is the energy density and p is the pressure. Raising one index, ^ has matrix form 


The continuity equation 


(-P 0 0 0\ 

0^00 
0 0^0 
Vo 0 0 y 


p' + 3-(p+4) =0 

a V & / 


(A31) 


(A32) 


may be obtained from the divergence-free condition on the stress tensor = 0. Letting w = pj(?p be the equation 

of state and dividing through by p, Eq. (IA32() is rewritten as 


-+3{l+w)— = 0. (A33) 

p a 

For the perfect fluids upon which most cosmological models are based, the equation of state is a constant, and we 
therefore find 


^ 

Next, there are two independent Einstein equations 



and 



(A34) 


(A35a) 


ca 


(A35b) 
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In the previous subsection on Cartesian coordinates we found the relations a = ta' and t = Utaoa~^. Given that 
changing to comoving spherical polar coordinates, or even physical spherical coordinates, leaves the time coordinate 
invariant we should expect the equation for t(r) to remain unchanged, as will be explicitly demonstrated below. Hence 
a' = in all coordinates under consideration, and Eq. (IA35al) becomes 


8 /ao\3(i+’") 

u'^c^aQ \ a J 


= 0 , 


(A36) 


which has solution 


a = ao 



2 

5+3u; 


UtCT 


(A37) 


Note that this solution looks rather different from the usual derivation of this sort, in particular, for a cosmological 
constant with w = —1 one normally obtains a = . The difference stems from the fact that we are parameterizing 

with respect to the affine parameter r rather than the coordinate time or, as is often done, the conformal time ry. 

The Cartesian parametrizations of null geodesics in FLRW is finally obtained by using the form of a(r) given by 
Eq. (IA37|) in Eqs. (IA24al) and (IA24bl) . Choosing the positive branch we find, for w — 1 


t = 


3(1 + w)c \3 




- 1/2 


/ 8 / 5 + 3r(; 

1 + \/ gTrpo 


3(l + '->) 

5+3™ 


UtCT 


- 1 


^0 


and 


X = 




- 1/2 


UfC / \1 + 3w/ \3 
with similar results for y and z. For ic = — 1, air) becomes 


,8 / 5 + 'iw 

1 + \/ giiPo 


1 + 3.1 
5 + 3.1 


UtCT 


- 1 


Xo, 


(A38) 


(A39) 


with 


and 


a = ag 






1 ± 


-Trpo UtCT 


+ to 



(A40) 


(A41) 


(A42) 


One of the fascinating things about FLRW spacetimes is that the peculiar motion (that portion of the motion that 
is not due to the Hubble flow) of timelike objects gets damped out, and eventually the physical motion is due solely 
to the Hubble flow. In other words, everything eventually becomes comoving. But the explicit parametric equation 
for the x-coordinate of null geodesics Eq. (IA42p shows that the same holds true for light in a universe dominated by 
a cosmological constant, w = —1. Asr^oo the null geodesic limits to comoving coordinate 


x(t —>■ oo) = Xo + 




(A43) 


The physical distance, of course, continues to increase such that a physical observer always measures c = 1 for the 
speed of light. Note that the scale factor when the light was emitted, ag, enters UtC through the null condition. This 
shows that in a cosmological constant dominated universe light is limited in the comoving distance it is able to travel. 
Thus all observers separated by a comoving distance 


B 




(A44) 


at r = 0 are not causally connected. 


































18 


4. Null geodesics in FLRW: comoving spherical coordinates 

Next, let’s consider spherical polar coordinates with line element given by Eq. (fTSl) . The null condition is 

- {rl + + sin^ = 0. (A45) 

and the geodesic equations are now 

^ + ’"c + sin^ = 0 (A46a) 


fc + 2—fc — rc{0^ + sin^ 0(p'^) = 0. 


(A46b) 


' + 2 \ — ] 6 — cos 9 sin 6ip‘‘ = 0 

rn a 


(A46c) 


Lp + 2\ — + - + cot 99 \ If = 0 

Tn a 


(A46d) 


From the null condition and a = a't, Eq. (IA46a|) leads immediately to 


i T — t — 0. 
a 


(A47) 


That the t equation is the same as that found previously is entirely expected; the change of coordinate system from 
Cartesian to spherical leaves the t coordinate unchanged, and hence the parametric equation for t should also remain 
unchanged. A first integral for ip can also be found immediately 


(f = UJo 


rcpao sin 9o 
Tra sin 9 


(A48) 


However, to find complete solutions it is advantageous to assume that 9 = 7r/2, in which case Eq. (jA46cl) also yields 
a first integral 


= 9n 


^Cq ^0 

Tea 


(A49) 


What this shows is that if we confine ourselves to the equatorial plane with 0o = 0; then we will remain in the 
equatorial plane. Finding the solution to the r equation Eq. (jA46b() is more difficult. In the limited case of purely 
radial null geodesics then 9 = ip = 0, and the equation reduces to 


Tc a 

+ 2 - = 0 , 

rr. a 


(A50) 


which is straightforward to solve. 

The task is somewhat less straightforward for non-radial null geodesics. As in the previous Minkowski example, we 
restrict to 9 = tt/2, whence the radial equation can be written 



2 4 4 


(A51) 


The complexity of the last term prevents us from finding an analytic solution in terms of Tc and a. Although Eq. (IA37I) 
provides an explicit expression for a(r), event this is too complicated to solve. On the other hand, this differential 
equation for is second order, but the null condition gives us a first order equation 


2 -2 I 

a r^ + 


2 4 4 


- 2=0 

a2 


(A52) 
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By substituting Eq. (IA37|) for a(T), one eventually finds the solution 


4t) = MtC" - ^ I 

Utcoofl + ow) \3 
This can also be written in terms of a(r) as 


1 ± 


8 75 + 3ri; 

o'^PO 


l + 3u 
5 + 3u 


UtCT 


-1 ± 


^0^ Cq 


aprcQ^o 


UtC 


rdr) = 


2 2 

UtC 


ao 


1 \ , aoTcgUr 


± 


1 + 3w \a(T)a(r) aoJ UtC^ 


UtC 


(A53) 


(A54) 


Knowing rc(T) explicitly, Eq. (IA48P can now be integrated to find an explicit parameterization of 
(/j(r) = tan“^ 


2 2 
c ui 


ao^’cn^o \ Utcao(l+3ui) ^3 




^ , /8 75 + 3u; 

1 ± \/ gTi-po 


1 + 3X. 
5 + 3x( 


utcr 


-1 ± 


aprcpUr^ 

uic^ 


=f tan 


Cco^O 


-ipo- (A55) 


Choosing an initial point (to,"Cco j ^/2, </2o) and initial velocity (uj,0, Wq) we then have an explicit parameterization 
for the curve. Furthermore, one can readily convert this to Cartesian coordinates with Xc(t) = rdr) cos((/j(t)) and 
yc(j) = r'c(T) sin(i^(T)). These become 


Xc{t) = 

Vc{t) = 


Utcao(l + 3w) \3 
2 

Utcao{i + 3w) \3 


Trpo 


, , 8 f5 + 3w 

1 ± \l 


l + 3u 
5 + 3u 


utCT] - 1 {ur^ cos ipo - rcpLOo sin (po) + Vcp cos po 




/8 / 5 + 8w 

1 ± \/ gTrpo 


1 + 3^ 
5 + 3x( 


UtCT ] - 1 I {ur^ sin Po + TcpCOo cos Po) + Tcp sin po- 

(A56) 


It can be useful to write down these main results in a more compact notation. We let 

2 


d = 


5 + 3w 


and A = Utc\ —npo- 


Then we have 


a = Go 




for the scale factor. Equations (IA41I) and (IA42I) become 


t = 


utd 


(1 - d)A 


l-d 




■ to 


and 


A 


(1 - 2d)A 

The spherical coordinate results Eqs. (IA54I) and (IA55I) become 


l-2d 




Xo- 


(A57) 


(A58) 


(A59) 


(A60) 


rc(r) = 


2 2 

UtC 


ao(l — 2d)A 




l-2d 


-1 ± 


^ 0 ^ Cq'^Vc 


f aprco^o 

V UtC 


(A61) 
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and 


(/j(r) = tan 


2 2 
C uf 


ao^’co‘^0 \ao{l -2d)A 
While Eqs. (IA56|) become 




l-2d 


-1 ± 


^ 0 ^ Cq'^Tc 

ulc^ 


=p tan 


Jco^O. 


+ (A62) 


Xc{t) = 
Vc{t) = 


ao(l — 2d)A 
d 

ao(l — 2d)A 






l-2d 


- 1 {Ur^ COS ^PQ - Tco Wo sin ipa) + cos pa 


l-2d 


- 1 {Ur^ sin(/9o + CcoWoCOStpo) + rco sinv^o- 


(A63) 


As discussed in the Minkowski example, if one is only interested in is the shape of the curve and doesn’t care 
whether the parameter is affine or not, then by rescaling with 


s = 


ao(l — 2d)A 


A 


l-2d 




(A64) 


we recover 


Xc(r) = s (ur, cospo - rcoC^osinipo) + rco cosipo 
ddr) = s [ur^ sin(/Jo + Tco^qCOSlpq) + Tco sint/^o, 

which is exactly the same as the spatial parametrizations in Minkowski spacetime. With an additional rescaling of 
t we have explicitly recovered the conformal flatness of FLRW in comoving spherical coordinates. This may seem 
like a very roundabout way to get to a result that could have been found more or less immediately from conformal 
arguments, however, the point of this approach is that we can find an explicit affine parametrization of the geodesics 
in these coordinates, and the fact that we explicitly recover the conformality of the solutions is an important check 
on the calculations. 


5. Null geodesics in FLRW: physical coordinates 

It is also advantageous to work in “physical coordinates” centered on the cloak, rather than co-moving coordinates. 
The reason for this is that for the cloak to remain of fixed size, the individual elements making up the cloak must 
overcome the Hubble flow by whatever internal forces that are keeping the cloak together. This means that although 
we may consider the center of the cloak to be comoving, other points on the cloak are not comoving. Therefore, 
a cloak of fixed physical size must actually be shrinking when described in comoving coordinates. Despite the fact 
that any coordinate system should be equally adequate for describing null geodesics, physical coordinates are may be 
more advantageous for describing a cloaking device of fixed physical size, and for describing the cloak parameters in 
a system in which all parts of the cloak are seen to be at rest with respect to the observer. 

Physical coordinates are reflective of the distances actually measured by an observer. A fixed comoving observer 
measures the physical distance from himself (at Vp = 0) to a galaxy at fixed comoving coordinate Vp to he d = a{t)rp. 
We let this increasing physical distance be the physical radial coordinate, so Vp = a{t)r, where r refers to the comoving 
coordinates. From this relation we have 


dvp = a'vcdt + advc = —rpdt + adrc, 


(A66) 


and therefore 


1 o! 

drc = -drp -yrpdt. 

a 


Plugging this into the FLRW line element Fq. there results 


ds^ = — 


a{t) 


2 P 


d^2 _ + drp + rpdi}^ 

aft) 


(A67) 


(A68) 
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Since the transformation from comoving to physical coordinates only involves transforming the r coordinate, it is 
expected that the equations for t and tp should remain essentially unchanged (we again assume equatorial motion). 
Calculating the geodesic equations confirms that expectation, with 


i = 


Uttto 

a 


if = LUq 



(A69) 


With motion constrained to the equatorial plane and these expressions for t and (p , the null condition becomes 




• 2 2 
a J'p 


2 4 


2 2 2 


= 0 . 


(A70) 


But it is straightforward to show that this same equation may be obtained by starting with the null condition Eq. (IA52I) 
in comoving coordinates and making the replacement ^ Tp/a. Therefore the solution may be found pretty much 

immediately from Eq. (IA53I) or Eq. (IA61I) . it is 


= ( 1 + dr 


\ 




(1 - 2d)A 


A 


l-2d 


-1 + 


Tpo (^i-p 4rp(,) \ Wq 


The equation for ip can now be directly integrated to find 

l-2d \ 


<p{t) = tan 


-1 


2 2 
ujc^ 


V (1 - 2d) A 


A 


-1 + 


fpo (Wr-p - Arp 
0,2^2 


ufc^ 


— tan 




- Ar„ 


Tpo^^O 


(A71) 


-po- (A72) 


Once again the Cartesian coordinates can be expressed as Xp = Vp cosp and yp = Vp sint/s. We find 


Xpir) 

yp{T) 


d 

(1 - 2d)A 



d 

(1 - 2d)A 




2d 



ij {{Urp 
1 ) (Kp 


/ y4 \ 

Arpg)cos(po -TpgUJo sin (fio) + ( 1 + ^po cospo 

Arpg ) sin po + ’’po'^o cos po) + + —t j Vp^ sin 


(A73) 

This is an interesting result that is substantially different to the usual behaviour an observer in Minkowski spacetime 
expects. As before, suppose that at r = 0 the curve makes its closest approach to Vp = 0, in which case the initial 
tangent to the curve has no radial component, Mpp = 0, and we can relate po to the angle pin at which the trajectory 
crosses the x-axis via po = (pm + "k 12. In both the Minkowski and comoving cases this allowed us to make the 
simplification 


x{s) = scospin — hsinpin 
y{s) = ssim^in + b cos pin, 


(A74) 


which carries with it the implication that if pin = 0 then the trajectory remains at fixed y and if pin = 7r/2 the 
trajectory remains at fixed x. But in the case of physical coordinates this is no longer the case. Indeed, for a trajectory 
with Ur = 0 and pin = 0, the y component continues to evolve as 


y(^) = ’’Po ( ^ + 7"^ 


2d 


I-2d 



(A75) 


Thus according to a physical observer, expansion of a congruence of null geodesics is inevitable and unavoidable. 


Appendix B: Coordinate transformations between comoving and physical coordinates 


The coordinate transformation between comoving coordinates (t,rc,0,p) and physical coordinates {t,rp,0,p) is 
given by 


(Bl) 


rp = a{t)r. 
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while the coordinate transformation from comoving to physical coordinates is given by the inverse of this, or 


rc = 




(B2) 


We may make manifest the invariance of a tensorial object under this kind of diffeomorphism of the manifold to itself, 
call it C : M ^ M and it’s inverse C~^. The associated differential dC : Tq^M —>■ Tq^M defines a map from the 
tangent space at a point Xc = qc in comoving coordinates to the tangent space at the related point Xp = qp = T{xc) 
in physical coordinates. The vectorial components at point Xc = q va comoving coordinates transform to the 
vectorial components Vq^ (where primed indices refer to physical coordinates) at point Xp = qp = C(p) in physical 
coordinates via the pushforward map given by the differential dC. This takes the form 


V 


v' 

Qp 


= Ju 


<lc 


(B3) 


where is the Jacobian matrix of C, evaluated at qc-, and is the matrix representation of the differential dC. Note 

that if Jp depends explicitly on the coordinate function Xc, then the components of the vector field v'^ will be 
expressed in terms of Xc, so we must furthermore use the coordinate transformation C to transform the coordinate 
function appearing in v'^ . 

Next consider the mapping of differential forms. A 1-form ujp in comoving coordinates must be pulled back to 
physical coordinates via the pullback map associated to C~^ rather than C itself. The covectorial components ujpiqc) 
at point qc in comoving coordinates transform to covectorial components uj'' (qp) (where primed indices refer to 
physical coordinates) at point Xp = qp = C{p) in physical coordinates via the pullback map given by the differential 
dC~^. This takes the form 


UJ,,>iqc) = L';^,\q^u}p{qc) (B4) 

where L^, is the Jacobian matrix of C~^ and is evaluated at the point qp. Since C~^ is the inverse of C, the two 
Jacobian matrices are related. In particular 


j/b-i 


K' = {j;) 


Xc—^C{Xc 


(B5) 


and the coordinate invariance of a tensorial object follows 

v(a;c) = 



d 

Qc 



= v„ 


‘‘^dx'i 

^C(x,) P 


7 = v{Xp). 


(B6) 


Now, if we want to transform the tensor components from comoving to physical coordinates, we apply 




yC{xc) 


(B7) 


A similar construction allows us to transform from physical coordinates back to comoving coordinates. The essential 
difference between a coordinate and corporeal transformation is that a coordinate transformation transforms both the 
components and the basis elements, as in Eq. (IB6I) . whereas a corporeal transformation acts only upon the components, 
while leaving the basis elements intact. 


Appendix C: Local frame transformations 

In curved spacetimes we often deal with coordinate patches where the coordinates are not orthonormal. This implies 
that the metric appears as something other than the Minkowski metric. It is desirable to be able to construct a locally 
flat frame, relative to which a vector, covector, or mixed tensor may be described. 

Begin with vectors. For a given choice of coordinates there exists a set of associated vector fields This set 

of basis vector fields span the tangent bundle over some portion of the manifold. Relative to this set of coordinate 
basis vector fields, any vector field may be decomposed as 


V = 


d 

dx°‘ ’ 


(Cl) 
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These basis vector fields are chosen out of convenience because they have a natural association to the chosen coor¬ 
dinates. However, even with the chosen coordinates it is possible to choose a different set of basis vector fields that 
will span the tangent bundle over the same portion of the manifold. Let us label this new set of basis vector fields 
{ba}, which is often called a tetrad or vierbein field. Since the basis sets {^f^} and {ca} both span the same vector 
(tangent bundle) space, we may express one set as a linear combination of the other. 


ga — Ga 


a 


d 


(C2) 


Generally, the set {ba} is not associated to a coordinate system in the same way that the set is. We therefore 
refer to the set {ba} as a non-coordinate set of basis vector fields. Since they both span the tangent bundle, the set 
{ba} is just as good a choice as {gf^} for decomposing a vector field. Thus the vector field V that was previously 
decomposed over the coordinate basis vector fields may now be decomposed over the non-coordinate vector fields 

V = l/^BA. (C3) 


One can think of the usual coordinate basis vector fields as a specially chosen tetrad field, subject to conditions 
discussed in section El 

Similarly, associated to the chosen coordinate system we have the set of basis covector fields (1-forms) {dx^}. This 
set of basis 1-forms spans the cotangent bundle over the same portion of the manifold as above. Relative to this basis, 
any 1-form lo may be decomposed as 


UJ = UJadx°‘. 


(C4) 


Again, these basis covector fields are chosen out of convenience because they have a natural association to the chosen 
coordinates. However, even with the chosen coordinates it is possible to choose a different set of basis covector fields 
that span the cotangent bundle over the same portion of the manifold. In particular, we will choose a set of covector 
fields {0^} that are dual to the vector fields chosen above such that O^es = Sg, just as for coordinate basis fields 
we have the dual relation dx^-^ = Sjj. Since the basis sets {0^} and {dx^} both span the same vector (cotangent 
bundle) space, we may express one set as a linear combination of the other. We let 

= e^^dx^. (C5) 

Since the set {0^} is not associated to a coordinate system in the same way that {dx^} is, we refer to {9^} as a 
non-coordinate set of basis covector fields. 

Because of the dual nature of both the coordinate and non-coordinate sets of basis fields, we have 


O^eg = {e^^dx>^)(eg^-^) = {e^^eg'^)dx>^^ = = e^^Cg^. 


dx° 


But since by duality 6^eg = 6g we have a condition on the coefficients, namely that 


A n xA 
2 „Bo'^ = Or. 


Now, since 


a 9 d d 


V = H^BA = F^eA“^ = (H^BA“) 


dx‘^ 


= 


dx° 


(C6) 


(C7) 


(C8) 


we have a way of transforming the coefficients of a vector field from the coordinate basis to the non-coordinate basis, 
and similarly for 1-forms. Thus we know how to transform both the coefficients and basis elements that describe 
vector fields and 1-forms. Transformation of mixed tensor fields follows straightforwardly in the usual manner. 

Given a line element 2-form ds^, it can be decomposed relative to either the coordinate basis fields, or any other 
set of basis fields. 


ds"^ = = Cji^eg^ gai3e^^dx°‘e^ pdx^ = gapdx'^dx^. 


(G9) 


In particular, we may choose a set of tetrad basis fields such that gAB is the Minkowski metric. However, that metric 
does not refer to Minkowski coordinates or even to the Minkowski coordinate basis covector fields. Instead, it refers 
to the non-coordinate basis covector fields {9^}. 

To summarize, we may decompose a vector field V, covector field us or mixed tensor % with respect to either 
coordinate or non-coordinate basis fields, and this can be done regardless of the choice of coordinates. A tetrad 
comprises a chosen set of basis vector fields at each coordinate point, but they are independent of the coordinate basis 
at that point. 
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Appendix D: Coordinate transformation between local frames 


Given Xai 3 ^''i^c) in comoving coordinates, we now understand this symbol to represent the coefhcients of the tensor 
object X decomposed relative to the sets of coordinate basis vector and covector fields (or some tensor product of 
them). From the previous discussion of local frames and tetrads, we can instead decompose x relative to a chosen 
frame where the metric is Minkowskian 


XAB^^ixc) = 


(Dl) 


On the other hand, in physical coordinates we have Xa'/ 3 '^ (^p) coefficients of the same object x decomposed 

relative to the coordinate basis vector and covector fields of physical coordinates. But we can just as well describe 
this relative to a locally flat set of basis vector fields as 


Xa'b 


C'D', 


a) — ^A' 


P'C’ 


P’ 


D' 


{Xp). 


(D2) 


We know how to transform the coordinate decomposition of x from comoving to physical coordinates, so the remaining 
question is how do we transform from the non-coordinate decomposition of x in a local frame of comoving coordinates 
to the non-coordinate decomposition of x in a local frame of physical coordinates? From Eqs. (iBTl) and (imi) this 
should be 


V (r \ — P P' pC pD' ja tP a B p V 1 T’^' 

XA'B' \^p) ® v' XaB y^cjOp, Op 

Xc^T{Xc) 

(D3) 

Combining the frame transformations with the coordinate transformations leads to 

Xa'B'^ ^ i^p) = Xab^^ 


(D4) 

3 

II 


(D5) 

and 

M^'c = 


(D6) 
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